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Abstract 

Using data from the FOCUS experiment we analyze the D + 7r~ and D°tt + in- 
variant mass distributions. We measure the D^f mass M D *o = (2464.5 ± 1.1 ± 
1.9) MeV/c 2 and width T D , = (38.7±5.3±2.9) MeV/c 2 , and the D* + mass M D „+ = 
(2467.6 ± 1.5 ± 0.76) MeV/c 2 and width r D „+ = (34.1 ± 6.5 ± 4.2) MeV/c 2 . We find 
evidence for broad structures over background in both the neutral and charged final 
state. If each is interpreted as evidence for a single L = 1, j q = 1/2 excited charm 
meson resonance, the masses and widths are M®, 2 = (2407 ± 21 ± 35)MeV/c 2 , 

T° 1/2 = (240 ± 55 ± 59)MeV/c 2 , and AfJ 2 = (2403 ± 14 ± 35)MeV/c 2 T+ /2 = 
(283 ± 24 ± 34) MeV/c 2 , respectively. 



Interest in charm spectroscopy has shifted from the ground states of (0~ and 
1~) cq mesons to the orbitally and radially excited states. In the limit of 
infinitely heavy quark mass, the heavy-light meson behaves analogously to 
the hydrogen atom, i.e., the heavier quark does not contribute to the orbital 
degrees of freedom (which are completely defined by the light quark). The 
angular momentum of the heavy quark is described by its spin Sq, and that 
of the light degrees of freedom are described by j q = s q + L, where s q is the 
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light quark spin and L is the orbital angular momentum of the light quark. 
The quantum numbers Sq and j q are individually conserved. The quantum 
numbers of the excited L = 1 states are formed by combining Sq and j q . For 
L = 1 we have j q = 1/2 and j q = 3/2. When combined with Sq they provide 
two j q = 1/2 (J=0,1) states, and two j q = 3/2 (J=l,2) states. In this paper 
these four states will be denoted by Dq, Di(j q = 1/2), D±(j q = 3/2) and 
D\. Heavy Quark Symmetry (HQS) predicts the spectrum of excited charmed 
states [1]- [5]. In the HQS limit, conservation of both parity and j q , requires 
that the strong decays D^j\j q = 3/2) — > D^n proceed only via a D-wave 
while the decays D^j\j q = 1/2) — > D^n proceed only via an S-wave. The 
states decaying to an S-wave are expected to be broad while those decaying 
in a D-wave are known to be narrow [6] [7]. Models predict that, when the 
heavy quark is the charmed quark, the physical states will have properties 
very close to those of the heavy quark limit. In the analysis described, we 
show the salient features of the D + tt~ and D°tt + invariant mass distributions 
and measure parameters of the well-established narrow states. We observe 
an excess of events in the mass interval 2.25GeV/c 2 to 2.40GeV/c 2 that is 
consistent with a broad resonance and must be included in the representation 
of the data to produce a good fit. 

The data for this paper were collected in the Wideband photoproduction 
experiment FOCUS during the Fermilab 1996-1997 fixed-target run. FO- 
CUS [8] [9] [10] is an upgraded version of experiment E687 [11] [12]. In FOCUS, 
a forward multi-particle spectrometer is used to investigate the interactions 
of high energy photons on a segmented BeO target. We obtain a sample in 
excess of 1 million fully reconstructed charm particles in three decay modes: 
D° — > K~n + , K~it + it + it~ and D + — > K~ir + n + . (The charge-conjugate states 
are implicitly included throughout the paper.) 

The FOCUS detector is a large aperture, fixed-target spectrometer with ex- 
cellent vertexing and particle identification. A photon beam, with an endpoint 
energy of ~ 300 GeV, is derived from the bremsstrahlung of secondary elec- 
trons and positrons. The charged particles which emerge from the target are 
tracked by two systems of silicon microvertex detectors. The upstream sys- 
tem [10], consisting of 4 planes (two views in 2 stations), is interleaved with 
the experimental targets, while the other system lies downstream of the tar- 
get and consists of twelve planes of microstrips arranged in three views. These 
detectors provide high resolution separation of primary (production) and sec- 
ondary (decay) vertices with an average proper time resolution of 30 fs 
for 2-track vertices. The momentum of a charged particle is determined with 
five stations of multiwire proportional chambers by measuring deflections in 
two analysis magnets of opposite polarity. Three multicell threshold Cerenkov 
counters [8] are used to discriminate between electrons, pions, kaons, and pro- 
tons. 
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1 Analysis Procedure and Results 



The L — 1 charm mesons were reconstructed via D + n~ and D°ix + combi- 
nations. The D° decays were reconstructed in the channels D° — > K~n + 
and D° — > i^~7r + 7r + 7r~. The D + decays were reconstructed in the channel 
D + — > A' _ 7r + 7r + . To obtain a clean sample of high statistics charm decays, 
the vertexing and particle identification cuts were optimized separately for 
each decay mode. The significance of separation between the production and 
decay vertex, Ijai (where t is the separation between the primary and sec- 
ondary vertex, and <jp is its error), was required to be greater than 5, 10, 
and 12 respectively for the three decay modes. The interaction vertex was 
formed from the D candidate, the bachelor pion and at least one additional 
charged track [11] and was required to be located within the target material. 
The pion and kaon candidates were required to have a Cerenkov identification 
consistent with the selected particle hypothesis. Further, we required that the 
decay D° — > K~7t + 7t + tt~ be reconstructed outside of target material and that 
| cos 6k\ < 0.7 for the D° — > K~n + decay, where 9k is defined as the angle 
between the D lab frame momentum and the kaon momentum in the D cen- 
ter of mass frame. Our starting samples for the decay modes with the above 
cuts are 210,000, 125,000 and 200,000 events respectively (see Figure 1 a-c). 
Combinations within ±2er of the nominal masses were retained as D candi- 
dates. Events with D° candidates coming from D* + decays were eliminated 
by applying a ±3a cut around the D* + — D° mass difference (see Figure 1 d). 

Figure 2a) shows the distribution of the invariant mass difference 

AM = M((K~7i + 7i + )n-) - M(K-7i + n + ) + M PDG (D + ) (1) 

where Mpdg(-D + ) is the world average D + mass [7]. Figure 2a) shows a pro- 
nounced, narrow peak near a mass M 2460 MeV/c 2 , which is consistent 
with the Df mass. The additional enhancement at M ~ 2300MeV/c 2 is con- 
sistent with feed-downs from the states D® and D^ decaying to D* + tt~ when 
the D* + subsequently decays to a D + and undetected neutrals. 

The mass difference 

AM+ = M((K-n + , K-n + n-n + )n + ) - M(K-n + , 7^+7^+) + 

+M PDG (D°) (2) 

spectrum (Figure 2b) shows similar structures to the AMq spectrum. The 
prominent peak is consistent with a of mass M ~ 2460MeV/c 2 . The 
additional enhancement at M ~ 2300MeV/c 2 is again consistent with feed- 
downs. 
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We fit the invariant mass difference histograms with terms for the DZ^, D% + 
peaks, Di and D\ feed-downs, combinatoric background and the possibility of 
a broad resonance. Fit terms were independent for each histogram except for 
specific systematic tests, and all fit parameters were allowed to float except in 
tests which are described below. 

The Dlf, D% + signals were represented with relativistic D-wave Breit-Wigner 
functions convoluted with a Gaussian resolution function (a = 7MeV/c 2 ). 
The a of the resolution function was determined by processing PYTHIA [13] 
events through the FOCUS detector simulation and reconstruction codes. 

The combinatoric background was represented by a continuum function dis- 
cussed below. The feed-downs were represented using line shapes determined 
by reconstructing simulated D*ir events as Dn. The masses and widths used 
for the D*n and Dn came from the PDG or from our fit to the D\ as de- 
scribed below. Only the amplitudes of the feed-downs were allowed to float in 
the fit. A relativistic S-wave Breit-Wigner function was used to represent a 
broad resonance contribution (motivated below). 

In order to determine functions for the combinatoric background, several stud- 
ies were performed. We studied the distribution of events in wrong sign combi- 
nations (the D* + (D°7t + )tt~ reflection from the D\ is very small), simulations 
where no L = 1 charm mesons are present, and data sidebands of the D + 
and D°. We found that in all these cases, the combinatoric background is well 
described by a single exponential beyond 2250 MeV/c 2 . Several functions with 
threshold characteristics (described in Section 2) were utilized to include in- 
formation below 2250 MeV/c 2 . Our final result is based on a function adapted 
from an E687 analysis [14] of excited D states 

exp{A + Bx){x-C) D (3) 

where x = AMo 5 +, and A, B, C and D are free parameters in the fit. (Care 
is required to limit the range of the C parameter so that the threshold term 
does not become imaginary.) With this function representing combinatoric 
background, we produced final results that were stable with consistently good 
confidence levels over a variety of fit ranges. No combinatoric shapes consistent 
with our background studies were able to describe either signal histogram 
unless we included a function representing a broad resonance. 

In order to illustrate the motivation for including the broad resonance, we 
show two representative fits performed without the broad resonance. The dis- 
tributions shown were fit with the D^ , parameters fixed to the world 
average values [7] (Figure 2 a,b), and with the D^ , D% + parameters allowed 
to freely float (Figure 2 c,d). Individual fit components, and an expanded view 
of the region around 2400MeV/c 2 are shown in the figure. In both cases, the 
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fit quality is unacceptable, even when the D\ parameters float to values far 
from the PDG values. Both fits indicate an excess of events between the 
signal and the feed-downs. We expect the background to be well described 
by a single exponential in this region, but the fit is unable to simultaneously 
describe the data at masses higher than the D\ peak and at masses lower than 
the D\ peak. Since the behavior of the combinatoric background is heavily in- 
fluenced by the events with invariant mass difference higher than the D\ peak, 
departures from the exponential form near 2400MeV/c 2 become evident. 

While we are unable to rule out the possibility that the excess is due to feed- 
down from higher mass charm states, we chose to describe the excess with an 
S-wave relativistic Breit-Wigner function centered roughly near the excess. 

In Figure 3 a-b) we show a fit to the data between 2120 MeV/c 2 and 5000 MeV/c 2 
that includes an S-wave relativistic Breit-Wigner in addition to previously de- 
scribed terms. Agreement is excellent with a fit confidence level of 22%. For 
self consistency, the D\ parameters measured in this fit are used to recompute 
the feed-down lineshape. When the histograms were refit using the new feed- 
down lineshape, the fit confidence level increased to 28% without a significant 
change in the returned fit parameters. The results of this last fit are shown in 
Table 1 together with PDG values where available. 

We find that the mass and width returned by the fit are increased compared 
to those reported by the PDG. Further, the yields and returned errors for the 
broad states indicate a significant excess is present. 

Although we are unable to distinguish between a broad state produced di- 
rectly via a Dq and the feed-down from a broad D\ state, we can make some 
qualitative comparisons. If the ratio of D\ to D\ production is the same for 
the charged and neutral modes, and the decays of these states are dominated 
by Dtt and/or D*7r, a meaningful comparison between the relative abundance 
of the feed-downs and the broad resonance to the D\ signal can be made. 
With these assumptions, one expects that the feed-down from the D\ and D\ 
narrow states to be larger relative to the D\ peak for D°tt + modes since the 
D*° has no D + channel. This is what we observe. We also find that the broad 
state contribution in the D°ti + mode relative to the peak is larger than 
the broad state contribution in the D + n~ relative to the D^ peak. This sug- 
gests some feed-down contribution to the broad state, perhaps from a broad 
D\ state (the search for a D*n broad resonance is being performed and will 
be included in a later publication on D*n states). 

Further, the fit parameters representing the broad S-wave state are statistically 
indistinguishable for both charged and neutral states. This is expected for 
broad states differing only by the flavor of the light quark and dominated by 
decay into a meson and a pion. 
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2 Systematic Checks 



Our systematic studies included a verification of the fit, fits using different 
functional forms for the background, different shapes for the feed-down, fits 
excluding the feed-down regions, fits over different regions of the data his- 
togram, a fit where we shifted our bin centers, a fit with the bin size reduced 
by a factor of 2, fits in which we excluded data where the background shape is 
expected to differ from that of Equation 3, and separate fits for particle and 
anti-particle distributions. All the contributions were added in quadrature (see 
Table 2) and are described in more detail below. 

The fitting algorithm was extensively tested by fluctuating the data histogram, 
comparing errors returned by the fit and the spread of parameters from re- 
peated trials. We have also performed repeated fits to histograms generated 
with the fit function. We observe that the goodness of fit is acceptable, that the 
central values are unbiased and that the errors correctly describe the variation 
of the central values over the trials. 

We split the sample into particle and anti-particle, producing two statistically 
separate data samples. These two samples were fit, and additional error (if 
any) was assessed until the parameters returned by the fit agreed with their 
average 

~^2(x — x avg ) = <7 stat + <7 extra (4) 

The deviations in the fit parameters returned by the tests described below were 
added in quadrature to the split sample estimate to assess a total systematic 
error. 

In addition to the modified E687 function, we fit the data with a pure expo- 
nential background function 

exp(A + Bx) (5) 

We also fit the data with a background function including a Gaussian term 
exp(A + Bx + Cx 2 ) (6) 

and we fit the data with a background function that was used by L3 [15] 

exp(A + Bx)/(l +exp(£> -x)/E) (7) 

We used feed-down functions based on PDG values for the D\ parameters. 
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Table 1 

Measured masses and widths for narrow and broad structures in D + ir~ and D°ir + 
invariant mass spectra. The first error listed is statistical and the second is system- 
atic. Units for the masses and widths are MeV/c 2 . 





^2 


D* + 

^2 


D*+ - Df 


D° 

U l/2 


D t/2 


Yield 


5776 ± 869 ± 696 


3474 ± 670 ± 656 




9810 ± 2657 


18754 ± 2189 


Mass 
PDG03 


2464.5 ± 1.1 ± 1.9 
2458.9 ± 2.0 


2467.6 ± 1.5 ± 0.76 
2459 ± 4 


3.1 ±1.9 ±0.9 
0±3.3 


2407 ± 21 ± 35 


2403 ± 14 ± 35 


Width 
PDG03 


38.7 ±5.3 ±2.9 
23 ± 5 


34.1 ±6.5 ±4.2 
25+ 8 7 




240 ± 55 ± 59 


283 ± 24 ± 34 



In addition we used feed-down functions based on our measured values for 
the D\ parameters. We also fit the entire histogram from 2030MeV/c 2 to 
5000 MeV/c 2 while excluding the feed-down region (2230 MeV/c 2 -2400 MeV/c 2 ) 
with both the E687 modified function and the L3 function, and we performed 
an additional fit with the E687 modified function where we exclude the region 
between 2120 MeV/c 2 and 2190 MeV/c 2 in addition to excluding the feed-down 
regions. 

We find that the data samples at very high (> 30) tjoi and high D momentum 
(P D > 70GeV/c) have a significantly different background distribution. We 
test the effect on our final result by removing these samples and by refitting. 

In order to determine the systematic uncertainty in our mass difference due to 
the mass scale of the FOCUS spectrometer, we measured the mass differences 
M(D*) - M(D) and M(i/j(2S)) - M(J/i/j). The quoted uncertainty is the 
additional contribution (added in quadrature to the statistical error) needed 
for our measurements to be in agreement with world average values. 

The contributions to the final systematic errors shown in Table 1 are listed 
in Table 2. The yields for both the narrow and broad states show a large 
variation depending on the fit considered. This is due to the wide range of 
background shapes investigated. Further, since the broad resonance is not 
fully contained in the fits, determination of the yield of the broad resonance 
depends on how much of the data histogram is included in the fit, and quoting 
a systematic error on this yield becomes problematic. Rather than quote a 
systematic error on the yield of the broad state, we looked at the statistical 
significance, Yield/5 (Yield), for each fit considered. In Figure 4 we show that 
the statistical signficance of our quoted result is a good representation of the 
fits tried. 
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Fig. 1. Invariant mass plots for: a)D° — > K ir + ; b)D° — > K 7r + 7r + 7r ; 
c)D + — ► K _ 7r + 7r + . Invariant mass difference plot for d) D* + — > D°7r + . 

Table 2 



Individual contributions to the systematic error. Units are MeV/c 2 





D*° 
^2 

Mass 


Df 
Width 


Mass 


D 2 *+ 
Width 


D° 

Mass 


D° 

Width 


^l + /2 
Mass 


D i,2 

Width 


D*+ - Df 

Mass 


l/o < 30 


0.160 


1.231 


0.134 


0.960 


0.926 


15.73 


0.050 


2.871 


0.294 


Part/Antipart 


1.67 





0.53 














31.4 





P D < 70GcV/c 


0.227 


0.705 


0.392 


1.983 


2.482 


8.509 


10.38 


2.500 


0.165 


Different Fits 


0.412 


0.272 


0.124 


0.693 


10.48 


43.95 


1.439 


8.635 


0.353 


Fit Regions 


0.376 


0.536 


0.174 


0.991 


1.571 


12.80 


1.209 


6.657 


0.315 


Feed-down tests 


0.633 


2.373 


0.262 


3.289 


32.71 


31.91 


32.45 


6.137 


0.443 


Binning tests 


0.442 


0.576 


0.113 


0.770 


6.584 


6.652 


6.380 


0.894 


0.550 


Mass Scale 


0.100 





0.100 





0.100 





0.100 





0.100 


Total syst. error 


1.94 


2.89 


0.76 


4.2 


35.1 


59.0 


34.7 


34.0 


0.91 
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Table 3 

Predicted mass differences with respect to the D meson compared to this result. 
The charged and neutral states are averaged. In the case of the broad state we 
compare our result to Dq only. Units are MeV/c 2 . 



Reference 


D 2 


n 

D\ 




D o 




jq = 3/2 


jq = 3/2 


jq = 1/2 


jq = 1/2 




3 P„ 


3 P, 


1 P, 
"1 


3 P„ 


This paper 


599 






538 


World Av. [7] 


593 


556 






Kalasnmkova ct al. (2002) [24J 


579 


562 


603 


564 


Di Pierro et al. (2001) [23] 


592 


549 


622 


509 


Ebert et al. (1998) [22] 


584 


539 


626 


563 


Isgur (1998) [21] 


594 


549 


719 


699 


Godfrey and Kokoski (1991) [3] 


620 


590 


580 


520 


Godfrey and Isgur (1985) [2] 


620 


610 


560 


520 


Eichten et al. (1980) [20] 


645 


637 


498 


489 


Barbieri et al. (1976) [19] 


428 


380 


339 


259 


De Rujula et al. (1976) [18] 


494 


464 


384 


374 



3 Conclusions 

FOCUS has measured the D + n~ and D°tt + mass spectra and provided new 
values for the masses and widths of the and Dp~ mesons (Table 1) with 
errors less than or equal to the errors on world averages. 

The D* 2 masses and widths measured are found to be higher than the world 
averages. We attribute the change to the inclusion of an underlying broad 
state. 

We find significant evidence for a broad excess which we parameterize with an 
S-wave resonance. Our results are consistent with a broad resonance occuring 
near 2400MeV/c 2 with a width of about 250MeV/c 2 in both the charged 
and neutral modes. We are unable to distinguish whether the broad excess is 
due to a state such as the Dq, predicted by HQS at M m 2400 and width 
~ 100 — 200MeV/c 2 , or due to feed-down from another broad state, such as 
the Di(j q = 1/2), or whether both states contribute. 

Evidence for L = 1 broad (S-wave) states has been previously presented in 
B decays by CLEO in the D*+n- final state [16], and BELLE [17] in the 
D* + 7i~, D + 7i~ final states. 

Our measurements are compared to theory predictions in Table 3. The D% 
masses are in good agreement with [21] and [23]. Reference [23], in addition, 
predicts a D^-Dq mass shift consistent with our evidence, while [21] predicts 
a shift with the opposite sign. 
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Fig. 2. The fit to the D + ir~ and D ir + mass spectra for the case where the D-wave 
mass and width are fixed to the PDG values, the background is described by Equa- 
tion 3, and no broad resonance is included, is shown in a) and b). The case where 
the D-wave mass and width are allowed to float in the fit is shown in c) and d). 
Note that the none of these fits provides a good description of the data between the 
feed-downs (~ 2300MeV/c 2 ) and the D* 2 peak (~ 2500MeV/c 2 ). In Figure 3, we 
show that the data are well described when a broad resonance is included in the fit. 
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AM GeV/c 2 AM + GeV/c 2 

Fig. 3. The fit to the D + ir~ and D°7r + mass spectra including a term for an S-wave 
resonance. The case with the mass and width for the Z?i(3/2) and D\ feed-downs 
fixed to the PDG values is shown in a) and b). The case with the mass and width 
for the D\{2>/2) feed-down fixed to the PDG values and for the D\ feed-down de- 
termined by fits in a) and b) is shown in c) and d). Notice the excellent agreement 
when the broad resonance is included (described in more detail in the text). 
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Fig. 4. Statistical significance of broad states signals in the D + ir (a) and D° 
(b) channels for various fits. 
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